Disentanglement in a quantum critical environment 
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We study the dynamical process of disentanglement of two qubits and two qutrits coupled to 
an Ising spin chain in a transverse field, which exhibits a quantum phase transition. We use the 
concurrence and negativity to quantify entanglement of two qubits and two qutrits, respectively. 
Explicit connections between the concurrence (negativity) and the decoherence factors are given for 
two initial states, the pure maximally entangled state and the mixed Werner state. We find that the 
concurrence and negativity decay exponentially with fourth power of time in the vicinity of critical 
point of the environmental system. 
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I. INTRODUCTION 

Entanglement is one of the most essential features in 
quantum mechanics [l[ and in recent decades has been 
focused by people in many fields of physics. Motivated by 
the progress of quantum information, entanglement has 
become a basic resource in the quantum technologies such 
as quantum teleportation and quantum cryptography (2j- 
4]. On the other hand, generally a realistic system is 
surrounded by an environment. The coupling between a 
quantum system and its environment leads to decoher- 
ence of the system. Thus, it is natural for us to consider 
the process of degradation of entanglement due to the 
decoherence. More recently, Yu and Eberly Q showed 
that two entangled qubits become completely disentan- 
gled in a finite time under the influence of pure vacuum 
noise. Surprisingly, they found that the behaviors of lo- 
cal decoherence is different from spontaneous disentan- 
glement. The decoherence effects take an infinite time 
evolution under the influence of vacuum while the entan- 
glement vanishes suddenly in a finite time. Some other 
researchers also investigated the process of disentangle- 
ment in the open quantum systems 6]- 1 9]. The problem 
of decoherence from spin environments was studied by 
Cucchietti et al while they considered the spin en- 
vironments consisting of N independent other than cor- 
related spins . 

In most of the previous studies, uncorrelated environ- 
ments are usually considered, and modelled by a reser- 
voir consists of harmonic oscillators. Although a collec- 
tion of harmonic oscillators is a well approximated mod- 
elling to represent the environment weakly coupled to 
system, however, in the practical situation, particles in 
the environment may have interactions with each other. 
Consequently, a problem comes out: How does the en- 
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tanglement evolves in a correlated environment? In this 
paper, we consider this problem and choose a correlated 
spin chain, the Ising model in a transverse field, as the 
surrounding system. Moreover, this surrounding system 
displays quantum phase transition (QPT) at some criti- 
cal point and thus it possesses the dynamic hypersensi- 
tivity with respect to the perturbation even induced by 
a single qubit jllj . 

As a quantum critical phenomenon, QPT happens at 
zero temperature, at which the thermal fluctuations van- 
ish. Thus, QPT is driven only by quantum fluctuation. 
Usually, at the critical point there exists degeneracy be- 
tween the energy levels of the systems when QPT hap- 
pens. Therefore, it can be excepted that, when we study 
the dynamic evolution of the system coupled to a envi- 
ronment with QPT, some special dynamic features will 
appear at the critical point. Quan et al [ll| have studied 
the decoherence induced by the correlated environment. 
It was shown that at the critical point of a QPT the de- 
coherence is enhanced. Following this work, Cucchietti 
et al [ijj discovered that the decoherence induced by the 
critical environment possesses some universality with the 
Boson-Hubbard model as an illustration. 

Now, we consider two spins coupled to the Ising spin 
chain in a transverse field, and the purpose is to reveal 
the effect of the correlated environment on the dynamic 
evolution of the two-spin entanglement. We will study 
different cases including two qubits and qutrits. More- 
over, we will consider cases that the two spins initially 
start from a pure maximally entangled state and a mixed 
Werner state [l3| . The 'sudden death' of entanglement 
is found to be a quite common phenomenon. 

This paper is organized as follows. In Sec. II, we intro- 
duce the model of two-spin system coupled to Ising spin 
chain with a transverse field. By exactly diagonalizing 
the Hamiltonian, we give expression of the time evolu- 
tion operator. In Sec. Ill, the analytical results of the 
concurrence [l4| of the two qubits are calculated to show 
the dynamics of entanglement. Numerical results are also 
given to illustrate the details of the dynamical behaviors 
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of entanglement. In Sec. IV, two qutrits are coupled to 
the Ising spin chain. The analytical and numerical re- 
sults of the negativity [HI, [l(| are given. At last we give 
the conclusion in Sec. V. 



From Eq. ((3|) and the units where % = 1, the time evolu- 
tion operator is obtained as: 



U(t) — e - i (-4 + 1 ) cr 0z t J~J gi^-ikx e ~itn k a kz e -i^-cr k:c _ ^ 
fe>0 



II. MODEL HAMILTONIAN AND EVOLUTION 
OPERATOR 

We choose the engineered environment system to be 
an Ising spin chain in a transverse field which displays a 
QPT. Two spins are transversely coupled to the chain. 
The corresponding Hamiltonian reads 
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where A characterizes the strength of the transverse field, 
g denotes the coupling strength between the Ising chain 
and the two spins, s\ and S2, erf (a = x,y, z) are the 
Pauli operators defined on the Z-th. site, and the total 
number of spins in the Ising chain is L = 2M + 1 . The 
Ising model is the simplest model which exhibits a QPT, 
and can be exactly calculated. 

In order to diagonalizc the Hamiltonian, firstly we no- 
tice that [s\ z + S2zi o~"] = 0, thus it is convenient to define 
an operator-valued parameter 



A = A+|(, l2 



S2z), 



(2) 



which is a conserved quantity. When we diagonalizc the 
Ising spin chain, the parameter A can be treated as a o 
numbcr with different values corresponding to the eigen- 
values of siz + S2z in the two-spin subspace. 

By combining Jordan- Wigner transformation and 
Fourier transformation to the momentum space [13] , the 
Hamiltonian can be written as [l8j 



k>0 



~ + A<T0z (3) 



where we have used the following pseudospin operators 
a ka (a = x,y,z) 

(T kx = dld i _ k + d- k d k ,(k = l,2,...M) 
a ky = -id' k d_ k + id-kdk, 



&kz = d\d k + d*_ k d- k - 1, 



002 



2dld Q - 1, 



(4) 



and d k ,d k {k = 0,1,2,...} denote the fermionic creation 
and annihilation operators in the momentum space, re- 
spectively. Here, 



A+2cos(27rfc/L) + 4 sin 2 {2nk/L), (5) 



i k = arcsm 



-2 sin 



'27rfc % 



fin 



(6) 



Having explicitly known the evolution operator, we now 
consider the entanglement dynamics of the two qubits 
and two qutrits. 



III. DYNAMICAL DISENTANGLEMENT OF 
TWO QUBITS 

A. The case with initial pure entangling state 

We investigate the dynamic evolution of two-qubit en- 
tanglement and assume that the two qubits initially start 
from a maximally entangled state. 
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(|00> + |11». 



(8) 



Here, |0) and |1) denote the spin up and down, respec- 
tively. The initial state of environment is assumed to 
be the vacuum state in the momentum space, namely, 
\iPe) = |0)fc=o <8>fc>o |0)fc|0)_fc, and the vacuum state 
|0)fe satisfies d k \0) k = 0. We may write a more general 
initial state of this composite system as 



|*(0)) = (o|00) + 6|ll))®|^) 



(9) 



From the evolution operator ([7]) , the state vector at time 
t is given by 

= a |00) ® U \tPe) + b |11) ® Ui He) , (10) 

where the unitary operator Uq and ~U\ can be obtained 
from the unitary operator U (t) by replacing operator A 
with number X+g/2 and \—g/2, respectively. 

Tracing out the environment, in the basis spanned by 
{|00) , |11) , |01) , |10)}, the reduced density matrix of the 
two-spin system is obtained as 



|a| 2 ab*F(t) 
p1 ' 2 ~ I a*bF*(t) \b\ 2 



^2x2, 



(11) 



where F(t) — {tPe\U\Uq \i/>e) is the decoherence factor, 
and ^2x2 denotes the 2x2 zero matrix. Now, the concur- 
rence [blj of the reduced density matrix can be readily 
given by 



C = 2\ab*F(t)\ = C Q \F(t)\, 



(12) 



where Co is the concurrence of the initial state. We see 
that the concurrence is proportional to the norm of the 
decoherence factor, and when the initial state is in a max- 
imally entangled state ijHJl, C — \F(t)\, namely, the con- 
currence is equal to the norm of the decoherence factor. 
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Let us consider the decoherence factor 

F(t) = (1> B \UlUo\il> B ) = Y[F k , 

k>0 



(13) 



where U n (n — 0,1) is generated from Hamiltonian H n 
with A = A„(a number). From the unitary operator ([7]) 
and the initial vacuum state, we obtain 



(0) 



k>0 



3 (!)1 2 



- cos(fi^t) sin(O^t) sin 6^] 

-sin 2 (^°^)sin 2 (^t)sin 2 (^-^)}^, 

(14) 

where fij^ and 6 k are obtained by replacing A with A„ 
in Eqs. §5§ and ©, respectively. Here, A = X+g/2 and 
Ai = X—g/2. This is one of our main results. We see 
that the zero mode (k = 0) has no contribution to the 
decoherence factor. Clearly, every factor F k is less than 
unit. So it can be well expected that in the large L limit, 
|-F(i)| will go to zero under some reasonable conditions. 

By carrying out similar analysis of Ref . [ll| , we intro- 
duce a cutoff number K c and define the partial product 
for the decoherence factor 



\F(t)\ c =l[F k >\F(t)\, 



k>Q 



from which the corresponding partial sum 



S(t)=]n\F(t)\ 



E 

fc>0 



InF, 



(15) 



(16) 



For the case of small k and large L, we have Q k 
1 2 — A„|, consequently 



(;n) 



sin 



16fc 2 7r 2 (Ap-Ai)" 
L 2 (2-A ) 2 (2-A!) 2 ' 



(17) 



As a result, if L is large enough and A — Ai is very small 
perturbation the approximation of S can be obtained as 

S(t) « ^2E(K c )(2-A Q y 2 (2-A 1 )- 2 

x{(A - Ai) 2 sin 2 (|2 - A | t) sin 2 (|2 - A x \t) 
+ [sin (|2 - A | t) cos (|2 - Ai| t) |2 - Ai | 
- sin (|2 - Ai| t) cos (|2 - A | t) |2 - A |] 2 }, 

(18) 

where 

E (K c ) = 4tt 2 K c (K c + 1) (2K C + 1) / (6L 2 ) . (19) 
In the derivation of the above equation, we have used 

n 

ln(l — x) w —x for small a; and ^ k 2 = n(n + l)(2n + 

fc=i 

l)/6. 
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FIG. 1: (a) Concurrence versus time t with different A in 
the case of weak coupling strength g = 0.1. The size of the 
environment is L = 300. (b) shows the cases of larger A. 




FIG. 2: Concurrence versus time with different environment 
size L = 200, 600 and 1000. The transverse field A = 4, and 
the coupling strength g = 0.1. 



For our two-qubit case, A = X+g/2, Ai = X—g/2. 
When A — > 2, and with a proper small g we have 



\F(t)l 



-if 



(20) 



with 7 = 2E(K C ) g 2 . Notice that \F(t)\ c is larger than 
\F(t)\ = C. Therefore, from the above heuristic analysis 
we may expect that when the parameter A is adjusted to 
the vicinity of the critical point A c = 2, the concurrence 
(or the decoherence factor) will exponentially decay with 
the fourth power of time. Moreover, for short times, from 
Eq. (fT4")l . the concurrence becomes 



C 



-rt 4 



(21) 



a(l) 
7 k 



)(nL 0) )W) 2 



with r = 1/2 Y, sin 2 

fc>0 

Now we resort to numerical analysis of the dynamical 
sensitivity and the concurrence decay. In the Fig. 1 (a) 
and (b) , we plot the concurrence versus time for different 
A. We find that in the vicinity of the critical point about 
Ae [2 — 0.3, 2 + 0. 3], concurrence decays monotonously 
with time. And extending the time range, however there 
are not the revivals of concurrence. Figure 1 (a) shows 
the cases of A < 2. We can see that concurrence for 
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FIG. 3: Concurrence versus time at the critical point A = 2 
with different coupling strength g. 



g = 25, concurrence will oscillate quickly and does not 
decay monotonously to zero any more. For the case of 
very large coupling g = 100, concurrence behaves as a 
weak oscillation near the initial value of C = 1. It can be 
expect that to the strong coupling limit of g, the concur- 
rence will stay at C = 1 without changing with time. The 
above behaviors remind us of the quantum Zeno effects 
in process of quantum measurement [21| . The phenom- 
ena shown in Fig. 3 is similar to the decay probability 
which can be suppressed by the increasing coupling be- 
tween system and measuring apparatus in quantum Zeno 
effects. 



B. The case of mixed state 



the case A = 2 decays more rapidly than other cases. It 
should be noted that, the dynamics of the two-qubit en- 
tanglement in Eq. (TT2"]) is absolutely determined by the 
decoherence factor in Eq. (fl4|) . thus from a theoretical 
point of view, the complete disentanglement cannot be 
realized in a finite time. When parameter A becomes 
larger than A c ,(g = 3,4 and 5), the numerical results of 
the concurrence are shown in Fig. 1 (b). The concur- 
rence oscillates with time, and collapses and revivals are 
observed. This is in contrast with the case of small A, 
where no revivals are found. 

The surrounding system displays a QPT near the crit- 
ical point, and there exists a competition between differ- 
ent order tendencies [l7l |. From another point of view, 
near the critical point quantum chaotic behaviors may 
emerge (l9| . For a system with quantum chaos, though 
it is prepared in identical initial state, two slightly differ- 
ent interactions can lead to two quite different quantum 
evolutions. In our system the decoherence factor can act 
as a fidelity and quantify the difference between the two 
states which are produced through two different evolu- 
tions. Decay of the fidelity can indicate the presence of 
the quantum chaos [2(J, and here the monotonous de- 
cay of the decoherence factor (concurrence) at the crit- 
ical point may be considered as a signature of quantum 
chaos. 

In Fig. 2, for weak coupling g = 0.1 and A = 4, the 
oscillation of concurrence is suppressed by enlarging the 
size of environment. The larger environment prevents 
the revival of entanglement. In the short-time region, we 
can see the larger size of environment will accelerate the 
monotonous decay of concurrence. From Eq. (|14[) . each 
factor Fk is smaller than 1 , thus it is reasonable that large 
size of environment will be more effective to suppress the 
factor F(t), and consequently suppress the concurrence. 

In Fig. 3, we consider the effects of coupling g on the 
dynamics of entanglement. At the critical point A = 2, 
we adjust g from a small one j = 0.1 to a strong one 
g = 100. It can be found that when we properly en- 
large the coupling, e.g. g = 1, the concurrence decays 
more sharply than the case g = 0.1. However, when we 
continue enlarging the coupling to about g > 10, e.g. 



Now, we study the dynamics of disentanglement of 
mixed entangled state and assume the two qubits being 
initially in a Werner state [l3j |. which is given by 

p s = P|$) <<I>| + I_£/ 4x4 , (22) 

where !<&) is the maximally entangled state given by Eq. 
©, the parameter P £ [0,1], and I^i denotes a 4 x 4 
identity matrix. This state is a mixed state except the 
extreme case of P = 1. Only when P > 1/3, the Werner 
state p s is entangled. 

We assume the initial state of the whole system p tot is 
in a direct product form as 

(hot = Ps <8 \^b) {iPe\ , (23) 

where \iPe) is the initial state of the environment. Af- 
ter the time evolution, we can obtain the reduce density 
matrix of the two-qubit system in the basis spanned by 
{|00) , |11) , |01) , |10)} as follows 

where the decoherence factor F(t) is the same as Eq. (|14|) . 
From Eq. ([2~4|) . the concurrence is derived as 

C = max|o,PMF| + -H. (25) 

When P = 1, it reduces to Eq. (p~2|) for the pure maxi- 
mally entangled state. While in the region 1/3 < P < 1, 
the concurrence vanishes when the decoherence factor 

1^1 < (P^ 1 - l)/2. (26) 

Thus there exists a finite disentanglement time t^, after 
which the entanglement is zero. According to the results 
of heuristic analysis in Eq. (f20|) . |P(t)| c ~ e~ 7 * , in the 
condition of weak coupling and A — > 2, we can approxi- 
mately give the disentanglement time 

>d °° (7 ln r^p)* ' < 27 » 



5 



O 



1 

0.8 
0.6 
0.4 
0.2 






--P=0.5 




— P=0.7 




— P=1 











4 6 
t 



8 10 



FIG. 4: Concurrence versus time at the critical point A = 2 
and coupling strength g = 0.1 for parameters P = 0.5, 0.7 
and 1. 



Then, the disentanglement time increases as the proba- 
bility P increases from 1/3 to 1. 

In Fig. 4, we also numerically calculate the concurrence 
versus time for different probabilities. For the mixed 
states corresponding to P = 0.5 and 0.7, disentangle- 
ment process takes only a finite time, while for the pure 
state case (P = 1), disentanglement is only completed 
asymptotically, and it will take an infinite time. Numeri- 
cal results are consistent with the above analytical results 
that the disentanglement time increases with the increase 
of P. 



IV. DYNAMICAL ENTANGLEMENT 
EVOLUTION OF TWO QUTRITS 

Now, we consider the case of two qutrits and use the 
negativity [i~5| to quantify entanglement. For the sys- 
tems with spin larger than 1/2, a non-entangled state 
has necessarily a positive partial transpose (PPT) ac- 
cording to the Peres-Horodecki criterion [15| . In the case 
of two spin halves, and the case of (1/2,1) mixed spins, a 
PPT is also sufficient. Vidal and Werner [l(| developed 
the Peres-Horodecki criterion and presented a measure of 
entanglement called negativity that can be computed ef- 
ficiently, and the negativity does not increase under local 
manipulations of the system. The negativity of a state p 
is defined as 



jv(p) = X>i, 



(28) 



where //, is the negative eigenvalue of p T2 , and T2 denotes 
the partial transpose with respect to the second subsys- 
tem. If M > 0, then the two-spin state is entangled. 
The negativity has been used to characterize the entan- 
glement in large spin system very well [22j]-[24|]. And 
by means of negativity, Derkacz et al. have studied the 
process of disentanglement in a pair of three-level atoms 
interacting with the vacuum Q. 



A. The case with initial pure state 

In a similar vein as the study of two-qubit case, we 
write a general initial state of the many-body system as 



|tf (0)) = (a |00) + b |11) + c|22» (8 \ipi 



(29) 



where |0), |1) , |2) denote the spin-one state with mag- 
netic quantum number 1, 0, -1 respectively. From the 
evolution operator (J7J), the state vector at time t is given 
by 

|¥(t)) = a |00) ® t/ +b |H) Ot/ilVs) 

+c|22) <g> U 2 \i(, E ) , (30) 

where the unitary operator Uq, J7i,and U 2 are obtained 
from the unitary operator U (t) by replacing operator A 
with number X+g, X and X — g, respectively. 

In the basis spanned by {|00), |11), \22), |01), 1 10), 
1 02), 1 20), 1 12), 1 21)}, the reduced density matrix of the 
two-qutrit system is 



H 2 ab*Fi(t) ac*F 2 {t) 
P1.2 = I a*bF?(t) \b\ 2 bc*F 3 (t) 
a*cF$(t) b*cF£(t) \c\ 2 

)Z 2x2 © ^2x2 © ^2x2, 



(31) 



where 



Fi(t) = WeWIUq \il> E ) , 
F 2 (t) = (1> E \U}U \il> E ) , 
F 3 (t) = (1>E\U}Ui\il> E ) 



(32) 



are the decoherence factors. 

The partial transpose with respect to the second sys- 
tem gives 



p{ 2 2 - diag(|or, |&r, |cH ©5iffiB 2 ffi B 3 , 
where the three 2x2 matrices 



(33) 



B 1 = 
B 2 = 
B. = 



aVFxit) 

a*bF?{t) 

ac*F 2 (t) 

a*cF*(t) 

bc*F 3 (t) 

b*cF£(t) 



(34) 



Then, from the above matrix p^ 2 2l one can obtain the 
negativity as 

N = \ab*F 1 (t)\ + \ac*F 2 (t)\ + \bc*F 3 (t)\. (35) 

For the maximally entangled state, a = b = c = 1/v3j 
and the negativity simplifies to 

A^=i(|F 1 (t)| + |F 2 (t)| + |F 3 (f)|). (36) 
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FIG. 5: (a) Negativity versus time with different cases of A = 
0.1, 1 and 2. The coupling g — 0.1 and the size of environment 
L = 300. (b) shows the cases of A = 3, 4 and 5. The highest 
one (solid line with up triangles) corresponds to the case A = 
5, and the lowest one (dashed line with points) corresponds 
to A = 3. 
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(a). When the coupling g is weak and A — > 2, the dynam- 
ical behaviors of the three decoherence factors in nega- 
tivity ([55]) are nearly identical. Each of the factors decay 
with time just as in Eq. ([2H)l . thus it can be understood 
that negativity also decays monotonously with time in 
the vicinity of A = 2. In Fig. 5 (b), we consider the cases 
of larger couplings. Comparing it with Fig. 1 (b), the 
behaviors of negativity have some differences with con- 
currence. More revivals are found in the behavior of the 
negativity, and they result from the linear superposition 
of the three decoherece factors. 

In Fig. 6, we numerically study the effects of different 
couplings g on the dynamics of negativity. Similar to the 
dynamic behaviors of the concurrence. With a properly 
large coupling such as g — 1, the decay of negativity 
will be much sharper. But very strong coupling (<? = 
15) will make negativity oscillate rapidly. To the strong 
coupling limit case of g = 100, negativity decays from 
the initial value JV = 1 to a steady value 1/3, which is 
different from the concurrence of the two qubits. Let us 
carry out the approximate analysis just like in the case 
of two qubits. We can obtain three partial sum Si, S 2 
and S3, corresponding to the three decoherence factors 
in Eq. (|32|) . which are similar to Eq. (p~8|) . When g — > 00 
and A — > 2, we have S 2 -> and Si = S 3 » -2E (K c ) t 2 
where E(K C ) is in Eq. (JTHJ) , thus negativity will decay 
sharply to a steady value of 1/3. We can see that different 
dynamic properties of the factors cause the behaviors of 
negativity shown in Fig. 6 is different from concurrence 
in Fig. 3. 



FIG. 6: Negativity versus time with different coupling 
strengths g = 0.1, 1, 15 and 100 at the critical point A c = 2. 

From the above equation, we can find the negativity is 
a linear combination of three decoherence factors. Also 
with the vacuum state of environment, the decoherence 
factors \F v (t)\ = {i/j E \UjUi \^ B ) are given by Eq.(H) 

by the replacements -> &f ~> ~> 

e k ) ^k ) -> 9 k ] - Hcrc ' F »( f ) denotes the three factors 
Fi(t), F 2 (t) and Fait). U]Ui correspond to U\U ,uIUq 
and U\U\ in the three factors Eq. (f3"2"|) . The parameters 
f4 n) and 9 { ™ ] {n = 0,1,2) can be obtained by substitut- 
ing Ao = X+g, Ai = A and A 2 = A— g into Eq. (JSJ) and 
©. 

During the similar analysis in the case of two qubits, 
we can also introduce the cutoff number K c and de- 
fine the partial product for the three decoherence fac- 
tors. Through the small k approximation, we can obtain 
the three partial sums corresponding to the three fac- 
tors. Therefore, under the condition of weak coupling g 
and A — > 2, in a finite time the three factors -Fi(t), -FbW 
and Fa(t) will decay exponentially with time in a similar 
form as Eq. ([20]), 

We numerically calculate the dynamics of negativity. 
In Fig. 5 (a), it shows the similar phenomena in Fig. 1 



B. The case of mixed state 

We then consider the mixed state, namely, the two- 
qutrit Werner state 

Ps = P\^)(^\ + l^I 9x9 , (37) 

where |$) is the maximally entangled state of two qutrits 
and |$) = (| 00) + |11) + |22)) /V3. Assume that the 
whole system is initially in p to t = Ps ® IV'e) {'4 ) e\- After 
time evolution operator in Eq. (JT]) , we can obtain the re- 
duce density matrix of the two qutrits at arbitrary time t. 
Then, we make the partial transpose with respect to the 
second system on the reduce density matrix, and obtain 

P T i% = idiag(l + 2P,l + 2P,l + 2P) 

®B X ®B 2 ®B 3 , (38) 
where the three 2x2 matrices 

B ^l{pWit) P ^ ) ) k={w} (39) 

From partially transposed reduced density matrix, the 
negativity is given by 

AA=i^ m ax{o,p(|F fe (t)| + i N )-H. (40) 
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Since |-Ffc(i)| < 1, the existence of nonzero negativity 
needs the parameter P satisfying the condition 1/4 < 
P < 1. From the above equation, we can also reads 
that the disentanglement occurs only when all the three 
factors satisfy \Fk(t)\ < (P~ x — l)/3. 

Furthermore, we study the case of a d-dimension 
Werner state being the initial state. Thus we give the 
initial state of the system as 

p v— -i i — p 

Ps = -J I") + —^p~ I d 2 xd^ (41) 

»,i=o 

where the basis vector \ii) is the eigenvector of s z = 
s iz + s 2z with the eigenvalue 2i + 1 — d. Then the initial 
state of the whole system is also performed by a direct 
product form as ptot = Ps <8> \^e) {iPe\ ■ After the simi- 
lar process mentioned in the former parts, we have the 
matrix p^ 2 2 denoting the reduce density matrix after the 
partial transpose over the second subsystem at time t, 
which is shown as: 

P 1 — P 

i,3=0 

= Idiag [1 + (d - 1)P, 1 + (d - l)P] dxd 

e . If ^ P M*)\ U2) 
ffi,<3 <lU^(t) ±f J' (42) 

where the decoherence factors Fi.j(t) — (ipE\UjlIi \"4>e) , 
and the corresponding time evolution operator Ui can be 
obtained from Eq. ((TJ) by replacing operator A with value 
X+g/2(2i + 1 — d), respectively. It is apparent that we 
should only focus on the 2x2 matrices and obtain the 
negativity as 

A r =igmax|o,P^ lJ (t)| + - i j , (43) 

from which we can see that negativity will be complete 
vanishes when all the norms satisfy |.fij(t)| < (P^ 1 — 
l)/d simultaneously. 

V. CONCLUSION 

In summary, we have studied the dynamics of entan- 
glement in a pure dephasing system. By making use of 



the concept of concurrence, we studied two qubits cou- 
pled to an Ising spin chain in a transverse field. When 
the two qubits initially started from a pure entangled 
state, we obtained the analytical results of concurrence 
which is just a simple product of the initial concurrence 
C(0) and the decoherence factor F(t). Thus the dynamic 
properties of concurrence is absolutely determined by the 
decoherence factor. Specially, in the case of weak cou- 
pling, the concurrence decays exponentially with time 
when A — ► A c . Moreover, we found the decay of deco- 
herence factor is of the form exp(— Tr 4 ), which is not a 
Gaussian form like in Ref. [ll[ and [T2] • Certainly this 
is due to the initial state of the environment we have 
chosen. 

Furthermore, when the two qubits arc initially in the 
Werner state, we have found that the complete disentan- 
glement takes place in a finite time just as the 'sudden 
death' of entanglement discovered in Ref. [f| . In [f| , due 
to the process of spontaneous emission, the sudden death 
of entanglement can occur in an arbitrary entangled state 
(pure or mixed). However, in our system with dephas- 
ing effects, when the two entangled qubits are in a pure 
state, there does not exist such a phenomena. 

We also considered two qutrits coupled to the Ising spin 
chain. When the qutrits initially start from a pure state, 
we have obtained the expression of negativity which is 
a linear combination of three decoherence factors. With 
weak coupling, negativity also decays monotonously in 
the condition A — > 2. When the qutrits are initially in 
a Werner state, the complete disentanglement could oc- 
cur in a finite time, and then the properties of negativity 
are the three decoherence factors. Indeed, the correlated 
environment, especially when QPT happens, greatly af- 
fects the decoherence and the disentanglement process. 
The entanglement decay in other environment which dis- 
plays a QPT [25| , or quantum chaos [2(| deserves further 
investigations. 
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